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ABSTRACT 


This  paper  describes  a  finite  element  approach  to  a  hybrid-mode  analysis  of  isolated 
and  coupled  microwave  transmission  lines.  Both  the  first  and  higher-order  ordinary 
elements,  as  well  as  singular  and  infinite  elements,  are  used  to  solve  for  the  eigenvalues 
and  eigenvectors  of  the  modes  propagating  along  the  line.  Once  the  field  distribution  in 
the  cross  section  of  the  line  is  known  the  characteristic  impedance  of  the  dominant  mode 
propagating  along  the  line  can  readily  be  obtained.  A  perturbational  approach  is  developed 
for  estimating  the  losses  due  to  conductor  and  dielectric  dissipation  and  computing  the 
attenuation  constant  for  each  mode. 

Lines  treatable  by  this  method  may  contain  an  arbitrary  number  of  arbitrarily 
shaped  conductors,  including  a  system  of  conductors  either  placed  above  a  single  ground 
plane  or  between  two  parallel  ground  planes,  and  inhomogeneous  dielectric  regions  that 
can  be  approximated  locally  by  a  number  of  homogeneous  subregions. 

The  results  obtained  using  the  finite  element  procedure  have  been  compared  for 
various  types  of  microwave  transmission  lines  and  have  been  found  to  agree  well  with 
available  data. 


I.  INTRODUCTION 


The  accurate  prediction  of  the  characteristic  impedance,  attenuation,  coupling, 
crosstalk,  etc.,  in  microstrips,  striplines  and  similar  transmission  lines  is  important  in 
microwave  and  millimeter  wave  integrated  circuits,  digital  circuit  design,  communication 
and  other  applications.  The  objective  of  this  paper  is  to  consider  microwave  transmission 
lines  with  rather  arbitrary  configuration  and  to  present  a  computer-aided  analysis  that 
allows  simple  and  accurate  calculation  of  its  design  parameters. 

In  the  past,  microwave  transmission  lines  have  been  thoroughly  investigated  by 
many  authors  and  many  approaches  to  analyzing  them  have  been  devised,  e.g.,  the  Green's 
function  techniques  [  1  ]-[€>],  conformal  mapping  [7]-[9],  variational  methods  [10],  [11], 
Fourier  transform  method  [12],  [li],  Fourier  integral  method  [14],  spectral-domain 
method  [  1 5]-[  17],  boundary  element  method  [18],  [19]  and  finite  element  method  [20].  All 
of  the  above  methods,  with  the  exception  of  the  last  two.  are  restricted  in  their  application 
to  transmission  lines  in  which  the  conductors  are  thin  strips  or  the  dielectric  inserts  have 
planar  interfaces  or  both  restrictions  apply.  In  contrast,  the  finite  element  method  (FEM) 
is  capable  of  handling  transmission  lines  with  rather  arbitrary  configurations,  since  the 
lines  treatable  by  this  method  may  contain  an  arbitrary  number  of  conductors  of 
arbitrary  shape  and  inhomogeneous  dielectric  regions  that  can  be  approximated  locally  by 
a  number  of  homogeneous  subregions.  Because  of  the  generality  of  the  FEM  approach  the 
finite  element  method  is  employed  in  this  papier  for  the  analysis  of  several  representative 
microwave  transmission  lines  of  practical  interest.  The  analysis  is  based  on  a  hybrid¬ 
mode  model  which  takes  into  account  coupling  between  TE-  and  TM-modes  propagating 
along  the  line. 


Both  the  first-  and  higher-order  ordinary  elements,  as  well  as  singular  and  infinite 
elements,  are  used  in  the  FEM  algorithm  to  solve  for  the  phase  constant  (eigenvalues)  and 
the  corresponding  field  distribution  (eigenvectors)  in  a  microwave  transmission  line.  The 


characteristic  impedance  of  the  dominant  mode  propagating  along  the  line  is  calculated 
from  the  so-called  power-current  definition  [2l]-[24].  A  perturbation  method  is  used  to 
compute  the  losses  for  each  mode  due  to  both  the  conductor  and  dielectric  dissipations. 

The  characteristic  impedance  and  loss  characteristics  have  been  calculated  for  some 
transmission  lines  of  interest  and  very  good  agreement  with  available  data  has  been 
obtained. 

II.  DESCRIPTION  OF  THE  METHOD 
1.  Formulation  of  the  Functional 

Consider  a  microwave  transmission  line  with  an  arbitrary  cross  section  consisting  of 
a  number  of  arbitrarily  shaped  conductors  and  inhomogeneous  dielectric  regions  which 
can  be  approximated  locally  by  homogeneous  subregions  (Fig.  1).  Let  us  assume  that  the 
line  is  uniform  along  the  longitudinal  (z)  axis  and  let  e  and  /jl  {~/j. 0)  denote  permittivity 
and  permeability,  respectively,  of  the  medium  in  each  homogeneous  subregion.  The 
electric  and  magnetic  fields  propagating  along  the  line  can  be  written  in  the  following 
form: 

E  =  ?  e~*:  .  R  —h  e~y'‘  (  1.1) 

where  y  -  j0  is  the  propagation  constant,  and  0  is  the  phase  constant. 

After  substituting  ( 1 )  into  first  and  second  Maxwell's  equation  we  obtain: 

V  x  /»  —  i  0(5  x  h  )  =  j  ueT  (1.2a) 

V  x  ?  -  /  0(:  x  e" )  =  —j  (A/jJi  (  1 .2b ) 

where  <x»=2v f  is  the  angular  frequency  of  the  wave  and  5  is  the  unit  vector  along  /- 
axis. 

Using  the  reaction  concept  [25]  we  derive  the  mixed-field  variational  formula  lor 
waveguide  phase  constants.  Multiplying  equation  (1.2a)  by  r  and  (1.2b)  by  h  and 


X 


Figure  I.  Cross  section  of  a  microwave  transmission  line. 


subtracting  I  1.2b)  from  (  1.2a)  we  obtain 


“  •  (  v  X/T  )  -  h  ■  (  v  x  c  )  —  f  ( h  x  c  )  j  &  —  :  ih  xc)-j&  =  j<dei'2  +  ijdi2)  (1.3) 
Using  the  relationship 

V'(3  xfl)  =  S-(Vx.4)-3-(yxfi)  (1.4) 

and  integrating  both  sides  of  (1.3)  over  the  cross  section.  S,  the  above  equation  becomes 

f  2 /  &:  ■  {c  x  h  )ds  —  f  V  •  (e  x  h  )ds  =  j  f  (  er2  +  fih  2)ds  ( 1.5) 

Js  s  s 

Finally,  af ter  applying  Stoke's  theorem  to  ( 1.5),  we  obtain 


a  J  (  €C2  +  fih2  )ds  —  x  h  )  •  ndl 

0=  — i - , - - - 

2  l :  •  [c  x  h  )ds 


(1.6) 


where  T  is  the  boundary  of  the  cross  section  S  and  n  is  the  unit  normal  at  the  boundary. 

Equation  (1.6)  is  the  variational  formula  for  0  which  is  stationary  for  small 
variations  in  electromagnetic  fields.  Reformulated  it  has  the  form: 

/*(  «’ :  +  jjJi  2 ) ds  -  —  f  f  •  (e  x  h  )ds  -  -L  6(c  x  h  )  •  ndl  =0  (1.7) 

Js  “  s  <*>  r 

In  the  general  case  electromagnetic  fields  ?.  /t  have  both  longitudinal  (z)  and 
transversal  ( 1 )  components 


0  =  c,  +  e.  f 


6  =  h ,  +  h.  z 


f  ransversal  fields  then  can  be  expressed  in  terms  of  longitudinal  components  e. ,  h. 


if  x  7  /?.  )  - 


where  V;  is  transversal  operator  and 


k  2  =  or  €fi  —  /32 

If  we  introduce  new  variables  [20] 

3=  ^  .  <t>  =  h:  .  «ir=  -v/ — 

0)  z  “  ft)  3 


relations  ( 1.9,10)  will  become 


t’r  =  1 


to  1 

Mr  —  7T 

<■  ** 


(f  xvi)-l  V,tf 

Mr 


A,  =  “  JT  |  V,  <t>  +  c>(2  X  V,tf)  ] 


Using  ( 1.12)  and  ( 1.13a,  b)  we  obtain 


e- 2  +  /jh  2  =  fj.0 


Hr<b2  + 


<i) 

c 


|Mr  Mr  +  32) l  V, d*l 2 


+  ^j32(^Mr  +32)l7fU»l2  +4erMr32f  X  V,  </>) 


2l(e  x  /T ) 

Ci) 


2j32Mr  I  V(  d>l  2  +  2B4€r  IV^I2 


+2/32  (/32  +  tMr  )5  (V,<&  X  Vr<*>) 


x  h  )  ■  n  =  — 


Mr(d>V,d>)  +  /32(t/>V,  t/t) 


n 


-j32  f  X  (iltV,<2>-<f>V,tf)  -n 

So  the  functional  (1.7)  finally  has  the  form 
k$  f  /xr<b2  +  Ji2\lr)ds  —f  r^r  IV ,  <l>\2ds  —f  r^.  j32!  V,»itl2ds 

*5  S  5 


( l.lb) 


—2  f  32rf  •  ( V(  t/t  x  Vt  <b)ds  — - 1 - -  Sie  x  h  )  •  ndl  =0  ( 1.17) 

\  t  IL  —  5Z)  T 


where 


r  =  (l  -32)/(^.Mr  ~  32) 


(1.18) 


If  the  last  term  is  omitted  from  the  functional  the  natural  boundary  conditions  of  the 
following  type  are  imposed  on  e  and  h  : 


1.  Tangential  electric  field  etang  is  equal  to  zero  on  the  surface  of  the  perfect 
electric  conductor  (PEC). 


2.  Tangential  magnetic  field  hlang  is  equal  to  zero  on  the  surface  of  a  perfect 
magnetic  conductor  (PMC). 


2.  Characteristic  Impedance  and  Dielectric  and  Conductor  Losses 

For  the  dominant  quasi-TEM  mode  the  characteristic  impedance  of  a  single  line  is 
calculated  using  the  so-called  power-current  definition  l2l]-[24] 


where  P  is  the  average  power  propagating  along  the  line 


8 


P  =  Re 

[  f  1 E  x  H*  Ls  ! 

=  /  Re 

1 e  x  h*  I  •  z 

[  s  '  •  1 

s 

(  J 

ds 


(2.2) 


and  I  is  the  current  induced  on  the  strip 


,mf'* 


xh,)dl 


(2.3) 


r5  is  the  contour  of  the  strip  and  n  unit  normal  on  the  surface  of  the  strip. 

Attenuation  constants  due  to  dielectric  and  conductor  dissipation  are  calculated  using 
a  perturbational  approach.  We  use 


<*d  = 


ac  — 


IP  '  2 P 

where  PD  is  the  power  dissipated  in  the  dielectric 

PD  =  tan  5  I  e  I  2ds  , 


(2.4) 


(2.5) 


tan  8  is  the  loss  tangent  and  Su  area  of  the  cross-section  covered  by  dielectric  and  e  is  the 
held  distribution  for  the  lossless  case. 


Pc  is  the  power  dissipated  in  the  conductors 


Pr  = 


hz2  +  I  n  x  h,  1 2 


dl 


(2.6) 


where  Rs  =  is  sur^ace  resistance  of  the  conductors  and  rr  is  the  boundary 


of  all  conductors  and  h, .  ht  are  held  distributions  for  the  lossless  case. 

Formulas  (2. 1-2.6)  can  be  expressed  in  terms  of  the  held  functions  <b  and  xb 
)2  . 

Mr  £  iVr<t>lzds  +  €r  B2  f  I  Vtrlf  I  2  ds 


P  = 


■Ti 


B  - 

€q  C 


1 


-  V"-.'-. 

i  *  •  ' 
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where  I",  is  the  path  from  the  middle  of  the  strip  toward  the  ground  plane. 


4.  Finite  Element  Solution 

According  to  the  FEM  procedure,  the  cross  section  of  the  line,  i.e.,  domain  S,  is 
subdivided  into  finite  elements  in  an  arbitrary  manner  provided  that  all  the  dielectric 
interfaces  coincide  with  the  element  sides.  Although  a  variety  of  different  elements  can 
be  chosen,  the  triangular  first-  or  higher-order  [26]  elements  are  adopted  in  this  study.  It 
has  been  shown  [20]  that  the  accuracy  and  efficiency  of  computation  can  be  substantially 
enhanced  via  the  use  of  the  high-order  elements  as  compared  to  the  case  where  only  the 
first-order  elements  are  employed,  if  the  same  number  of  nodes  is  used.  On  the  other 
hand,  the  use  of  a  large  number  of  simple  elements  is  clearly  advantageous  where  a 
complicated  boundary  shape  needs  to  be  modeled. 

For  ordinary  triangular  elements  the  area  coordinates  (Fig.  2)  are  defined  as  follows: 


&  = 


1 

2  A. 


(p,  +  </,  x  +  r,  y  )  ,  z  =1.2.3 


(4.1) 


Pt 


x'+\ 

y«+i 

xt  +2 

y<+2 

=  y,  +2  -  y,  +i  =  -{x<  +2  -  •*«  +i ) 


(4.2) 


where  \  is  the  triangle  area.  Coordinate  lines  =*  const,  are  shown  in  Fig.  3. 

For  the  high-order  elements  the  nodes  are  chosen  as  a  set  of  regularly  spaced  points 
(  Fig.  4).  For  the  element  of  order  M  the  number  of  nodes  is 


Aft'  =  (M  ■  (A/  + 1)3/2 


(4.3) 


which  is  equal  to  the  number  of  terms  in  a  complete  polynomial  in  x.y  of  order  M. 

Within  each  triangle  unknown  functions  6  and  xl)  are  expanded  in  terms  of  trial 

f unctions,  «,  and  the  nodal  potentials,  <1>; ,  'F, ,  i  -  1 . MV,  where  MV  is  the  number  of 

nodes  per  element,  as  follows: 


Figure  2.  (a)  Any  interior  point  P  defines  a  splitting  of  the  triangle 
into  subtriangles,  (b)  Area  coordinates  measure  relative  distance 
towards  each  vertex  from  the  opposite  side. 
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f  igure  4.  Point  placement  in  a  triangle  for  M  -  4.  showing  the  triple  and  single  index 
numbering  schemes  as  well  as  the  numbering  of  triangle  vertices. 
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For  an  ordinary  element  of  order  M.  the  functions  a,  are  complete  polynomials  in 
x,y  of  order  M  and  satisfy  the  following  conditions 

a,  —  1  for  x  =  X[  y  =yt  ( node  l  ) 

al  =  0  for  x  =  y j  y  =  y;  _  j  /  (4.5) 

In  accordance  with  the  triple  index  ordering  shown  in  Fig.  4,  trial  functions  can  be 
expressed  in  the  following  form 

cw,  =  ^  (*,)/»«  WP«  U.6) 

where 
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After  introducing  the  expansion  (4.4)  the  field  components  become 
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For  the  functional  (1.17)  we  need  to  calculate  the  following  terms: 
/'Vf  <t  I  “  ds  ,  I  V;  \ltl2  ds  ,  £  ~  *  ( V,  tit  x  Vi  <t>)  ds ,  and  ^  <b2  ds  ,  £  \lr  ds  where  Ae 

is  the  surface  of  the  element.  After  some  manipulation  they  can  be  expressed  in  the  form 
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where  S.  T  and  U  are  the  finite  element  matrices.  The  above  matrices  are  given  by 


(5  ]  =  [Q  i  ]  cot  6t  +  [Q ,]  cot  Bz  +  [Q  3]  cot  9j 


(4.14) 
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cot  9t ,  i  —  1,2,3  are  the  included  angles  and  matrices  Qt .  i  —  1,2,3  are  symmetric  and  do 
not  depend  on  triangle  shape.  The  elements  of  the  matrix  U  that  takes  into  account 
coupling  between  TE  and  TM  modes  are 

r -mn  1  V*  /*  (fam  3 an  (fan  ,  m  )• 

L  ds •  i.~Mv  «“-> 

This  matrix  is  antisymmetric  and  does  not  depend  on  triangle  shape. 

The  elements  of  matrix  T  have  the  form 
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(4.17) 


This  matrix  is  also  symmetric.  The  functional  ( 1.17)  now  has  the  form 

*<T  |  a,  <J>r  T<t>  +  JZ2Vr  TV  J  -  Tfxr  Or  5<t>  -  t%.  BzVt  S  V  ~2B2rVTU<&  =  0(4.18) 

where  <t>.^  are  column  matrices  of  nodal  potentials  and  <J>r ,  ^  are  their  transpose 

values. 

Unknown  nodal  potentials  are  found  by  differentiating  the  functional  (4.18).  The 
following  system  of  FEM  equations  is  then  obtained 
.v/v  /  ,  A/V 
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=i  1  •  i=i 

/  =  1.2.  •  •  •  A/V 

A/V  |  |  VA' 
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which  can  be  rewritten  in  the  form 
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where  .V  is  the  column  matrix  of  all  nodal  potentials  (¥, ,  <!>,  ),  A  =  J.  Equation 
(4.20)  is  an  ordinary  matrix  eigenvalue  problem  whose  eigenvalue  and  eigenvectors  are  k2 

and  ^  ,  respectively.  The  normalized  frequency  ^  is  related  to  k0 2  in  the  following 


(4.21) 


We  can  use  the  FE-matrices  to  express  the  other  quantities  of  interest  for  each 


element. 


1.  Power  propagating  along  the  line: 
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2.  Power  dissipated  in  dielectrics 
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3.  Electrostatic  energy  per  unit  length 
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4.  Current  induced  in  side  1-2  of  the  triangle 
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where  QIX.QI  i.  Q/ 3  are  the  column  matrices  of  dimension  MV  and  have  the  following 
form 
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5.  Voltage  induced  in  side  1-2  of  the  triangle 
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b.  Conductor  losses  in  side  1-2  of  the  triangle 
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+  py  cot  e,  cot  e2  +  tr2vr  py  cot  :0; 

where  ^,2  is  The  length  of  the  side  I -2  and  P0  —  Pb  are  the  FE-matrices  which  do  not 
depend  on  the  triangle  shape 
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5.  Special  FEM  Elements 

In  order  to  enhance  the  accuracy  and  efficiency  of  the  presented  FEM  procedure 
special  elements  as  singular  and  infinite  elements  have  been  used. 

When  there  are  field  singularities  caused  by  edges  in  the  cross  section  of  the  line,  a 
very  fine  mesh  of  first-  or  high-order  elements  is  required  to  obtain  an  accurate  solution. 
However,  in  order  to  improve  accuracy  and  reduce  the  number  of  necessary  nodes, 
singular  elements  [27],  [28]  are  used  in  this  study.  A  singular  element  with  associated 
triangular  polar  coordinate  system  (p.  a  )  is  shown  in  Figure  5. 


The  global  Cartesian  coordinates  are  then 


.v  =  x ,  +  p  [(.v,  —  x  , )  +  a  (,t3  —  ,r2)]  (5.1 ) 

y  =  v ,  +  p  [(y  2  -  v , )  +  a  ( v  3  —  y  2 )] 


The  potential  distribution  function.  <2>  or  tit,  that  takes  into  account  the  held 
singularity  at  the  node  1  (p  =  0)  has  the  following  form: 

<*>  =  I  *,(P)  t  (5.2) 

i*0  k =0 

where  R,  (p)  is  a  polynominal  of  the  form  {28] 

R,  (p)  =  ato  +  ai  \PK  +  atz  PX+1  "*■  +au  (5.3) 


and  A  is  a  coefficient  chosen  in  accordance  with  the  Meixner  edge  condition  [29].  L*(  a )  are 
Lagrange  interpolation  polynomials: 
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are  nodal  potentials  for  the  nodes  regularly  spaced  in  the  triangle  (Fig.  6). 
elements  are  compatible  with  high-order  ordinary  elements. 


(5.4) 

These 


Although  open-type  transmission  lines  can  be  treated  by  the  conventional  finite 
element  method  with  the  shielding  far  away  from  the  region  of  interest,  the  computation 
efficiency  can  be  substantially  improved  if  infinite  elements  [30]  are  used.  Consider  first 
the  line  consisting  of  a  number  of  conductors  arbitrarily  placed  between  two  parallel 
ground  planes  ( Fig.  7a). 

I  he  entire  domain  can  be  divided  into  the  near-held  region  ( n.l  . ).  which  is  the  region 
of  interest,  and  the  far-held  region  (f.f.).  which  is  unbounded.  They  have  a  common 
boundary  referred  to  as  the  far-held  boundary  (f.f.b).  The  near  held  (n.f.)  region  is  then 
divided  into  finite  triangular  elements  in  the  usual  manner,  while  the  f.f.  region  is  divided 
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I 

I 

into  infinite  elements.  Each  infinite  element  has  one  common  side  (lying  on  f.f.b)  with  an  j 

I 

ordinary  element  and  two  sides  parallel  to  the  x-axis.  | 

Consider  the  infinite  element  with  nodes  1  and  2  (Fig.  7b)  and  introduce  normalized 
coordinates 

£  =  — .  T)  —  — — — .  x,  <  x  <  oo,  y ,  ^  y  <  y,  (5.5) 

x\  y 2~y  i 

Thus,  in  the  normalized  coordinates  rj  the  potential  function  <t>  can  be  written  in  terms 
of  the  nodal  potentials  as: 

£  /,"(*)<&,  (5.6) 

i=0 

where  L-u  are  Lagrange  interpolation  polynomials  and  <&,  are  corresponding  nodal 
potentials. 

The  next  type  of  open  region  problem  is  a  system  of  conductors  above  a  ground 
plane.  The  procedure  is  similar  to  that  in  the  previous  case.  The  entire  domain  is  divided 
into  n.f.  and  f.f.  regions  (Fig.  8a).  The  n.f.  region  is  divided  into  the  usual  triangular 
mesh  and  the  f.f.  region  is  divided  into  infinite  elements  of  two  different  types.  1  and  II. 

The  infinite  element  type  I  is  the  element  previously  discussed.  Consider  now  the  infinite 
element  of  type  II  with  nodes  1,  2  and  radial  sides  intersecting  at  point  ( x0  y0  )  (Fig.  8b) 
and  introduce  triangular  polar  coordinates  p.o  which  are  related  to  the  global  Cartesian 
coordinates  by  the  relations 

*  =  x0  +  p  [  (x  x—x0  ,  +  a(xz— x  , )  ]  (5.7) 

y  =  v0  +  p  [  (v  ,-v0 )  +  o( v j  v , )  ] 

In  the  similar  way  as  for  the  element  type  I  the  potential  distribution  within  element 


type  II  is  approximated  by 
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Figure  8a.  I  me  consisting  of  a  system  of  conductors 
placed  above  a  ground  plane. 
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Both  infinite  elements  type  I  and  type  II  are  compatible  with  high-order  ordinary 
elements  of  order  M. 

Coefficient  a  is  varied  as  well  as  the  position  of  the  f.f.b  until  the  convergence  of  the 
solution  is  obtained. 


III.  NUMERICAL  RESULTS 

On  the  basis  of  described  finite  element  procedures,  two  different  computer  software 
packages  have  been  developed  for  solving  the  problem  at  hand.  The  first  of  these  is 
employed  for  semiautomatic  mesh  generation  of  first-  or  high-order  ordinary  elements 
including  singular  and  infinite  elements  if  necessary.  The  second  package  solves  for  the 
dispersion  characteristics  and  the  field  distributions  of  the  propagating  modes.  Using  the 
obtained  field  distribution  for  each  mode  it  then  calculates  the  characteristic  impedance  of 
the  line  and  attenuation  due  to  conductor  and  dielectric  losses. 

Different  types  of  microstrip  lines  have  been  analyzed  and  some  of  the  results  are 
presented  in  this  study. 

In  Fig.  9,  the  generic  cross-section  of  the  shielded  microstrip  line  is  presented.  The 
width  of  the  line  is  2w  and  the  thickness  t.  Height  of  the  substrate  is  h  and  the  height  of 
the  shielding  is  b.  The  width  of  the  shielding  is  2a. 

For  the  special  case  when  the  permittivity  of  the  substrate  is  «,  =  4^,  dispersion 
characteristics  of  the  first  four  propagating  modes  are  presented  in  Fig.  10.  Values  of  /32 
are  plotted  versus  normalized  frequency  (oi h  /c  )2.  It  should  be  pointed  out  that  for  the 
thickness  of  the  substrate  h  -  1mm  actual  frequency  for  (to/i  c  )2  =  1  is  about  50  GHz. 
The  dominant  mode  is  a  quasi-TEM  mode.  Another  mode  which  propagates  from  the 
zero  frequency  is  the  surface-wave  mode.  The  other  two  modes  are  waveguide  modes.  In 
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Fig.  1 1.  the  results  for  the  characteristic  impedance  of  the  dominant  mode  are  presented 
up  to  a  frequency  of  40  GHz.  In  this  case  permittivity  is  c  =  4^  and  the  height  of  the 
substrate  h  =*  0.5mm. 

For  the  same  line  the  results  for  the  attenuation  constant  due  to  conductor  losses  are 
presented  in  Fig.  12  and  for  the  attenuation  constant  due  to  dielectric  losses  in  Fig.  13. 

Also,  the  coupled  lines  have  been  analyzed.  The  generic  cross-section  of  one  half  of 
the  line  is  presented  in  Fig.  14. 

Both  even-  and  odd-modes  have  been  analyzed.  In  Fig.  15.  the  dispersion 
characteristics  of  the  first  four  propagating  even-modes  are  presented  while  in  Fig.  lb  the 
odd-modes  are  plotted. 


IV.  CONCLUSIONS 

In  this  study  a  full-wave  analysis  of  microwave  transmission  lines  has  been 
performed.  The  analysis  is  based  on  a  hybrid-mode  model  which  takes  into  account  the 
coupling  between  TE-  and  TM  modes  propagating  along  the  line.  The  finite-element 
method  has  been  used  to  solve  for  the  propagation  constant  (eigenvalues)  and  the  field 
distribution  (eigenvectors)  of  each  mode.  High-order  ordinary  elements  as  well  as 
singular  and  infinite  elements  have  been  employed.  Using  the  obtained  field  distribution 
the  characteristic  impedance  of  the  dominant,  quasi-TEM,  mode  has  been  calculated. 
Finally,  a  perturbational  approach  has  been  used  to  compute  the  attenuation  constant  due 
to  conductor  and  dielectric  losses. 

A  general  computer  software  package  has  been  developed  for  semi-automatic  mesh 
generation  and  for  computation  of  the  propagation  and  attenuation  constants  and  the 
characteristic  impedance.  Some  of  the  results  for  isolated  and  coupled  microwave 
transmission  lines  are  presented  and  a  good  agreement  with  available  data  has  been  found. 
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Figure  11.  Characteristic  impedance  of  shielded  micro-strip  line 
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Figure  14.  Generic  cross  section  of  one  half  of  shielded  coupled  microstnp  line. 
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